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FOREWORD 


This  report  summarizes  the  results  of  research  performed  by  the  Aeronautical  Systems 
Division,  Wright-Patterson  AFB,  OH.  The  work  was  performed  between  November  1976  and 
March  1978,  and  the  Project  Engineer  for  this  effort  was  G.  B.  Bennett. 

The  work  was  sponsored  by  the  JTCG/AS  as  part  of  the  3-year  TEAS  (Test  and 
Evaluation  Aircraft  Survivability)  program.  The  TEAS  program  was  funded  by  DDR&E/ 
ODDT&E.  The  effort  was  conducted  under  the  direction  of  the  JTCG/AS  Vulnerability 
Assessment  Subgroup,  as  part  of  JTCG/AS  Work  Unit  VA-6-02F,  Development  of  Air- 
craft Preliminary  Design  Assessment  Metfwdology  '.  — 

This  report  presents  and  summarizes  two  sorting  algorithms,  PRIME  Sort  and  PDQ 
Sort,  developed  in  the  Aeronautical  Systems  Division  Computer  Science  Center  in  support 
of  vulnerability  assessment  computer  programs  in  use  by  the  Deputy  for  Development 
Planning. 


NOTE 

This  technical  report  was  prepared  by  the  Vulnerability  Assessment  Subgroup  of 
the  Joint  Technical  Coordinating  Group  on  Aircraft  Survivability  in  the  Joint 
Logistics  Commanders'  organization.  Because  the  Services’  aircraft  survivability 
development  programs  are  dynamic  and  changing,  this  report  represents  the  best 
data  available  to  the  subgroup  at  this  time.  It  has  been  coordinated  and  approved 
at  the  JTCG  subgroup  level.  The  purpose  of  the  report  is  to  exchange  data  on  all 
aircraft  survivability  programs,  thereby  promoting  interservice  awareness  of  the 
DoD  aircraft  survivability  program  under  the  cognizance  of  the  Joint  Logistics 
Commanders.  By  careful  analysis  of  the  data  in  this  report,  personnel  with  ex- 
pertise in  the  aircraft  survivability  area  should  be  better  able  to  determine 
technical  voids  and  areas  of  potential  duplication  or  proliferation. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  this  PAGE  f'FTi.n  Cm  Colored) 

( REPORT  DOCUMENTATION  PAGE 


| JTCG/ ASf78-V-004 


TITLE  (end  Submit)  i 

PRIME  & PDQ  Sorts  - Efficient  Minimal  Storage 
Sorting  Algorithms 


7.  AUTjIOJUaL--  

)>'  oy  R.jHilbrand 

*.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Aeronautical  Systems  Division  , , j 

ASD/ADDP  0 " * 

Wright-Patterson  AFB,  Ohio  45433 

<1.  CONTROLLING  OFFICE  NAME  AND  AODRESS 

JTCG/AS  Central  Office 

Naval  Air  Systems  Command,  AIR-5 204 J 

Washington,  D.  C.  20361 


READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

NO.  S'  RECIPIENT'S  CATALOG  NUMBER 


Final  ^ •£ 

NovaPBWR7 


SPt» 

76-Mard 


. * $ 


CONTRACT  OR  GRANT  NUMBERS) 


program  element,  project,  task 

AREA  A WORK  UNIT  NUMBERS 


Work  Unit  VA-6-02F 


) Aug*PP79 

fe  NUMBER  OP  S 

54 


MONITORING  AGENCY  NAME  ft  ADDRESSfff  dllisrsnt  from  Controlling  OHIc»)  I 15  SECURITY  CLASS,  (oi  rhfa  report; 


\i)  L i f 


Unclassified 


tSa.  OECLASSIPICATION/ DOWNGRADING 
SCHEDULE 


[IS.  DISTRIBUTION  STATEMENT  (oi  this  Rsport) 


Approved  for  public  release;  distribution  unlimited.  Statement  applied  August  1979. 


17.  DISTRIBUTION  STATEMENT  (o I Ml.  ab.lr.el  entered  In  Block  10.  II  dllloront  from  Report)  

/ (0 J PRIMS  and  PDQ  Sorts  - Efficient 

Minimal  Storage  Sorting  Algorithms. 

IS.  SUPPLEMENTARY  NOTES 


IS.  KEY  WOROS  (Continue  on  rororoo  tide  II  nocoooory  end  Identity  by  block  number) 

Minimal  Storage  Sorts  Sorting  algorithms 

Computer  programs 
Vulnerability  Assessment 
Missile  endgame  models 

SO.  ABSTRACT  (Continue  on  rover  re  tide  It  nocoooory  end  Identity  by  block  number) 


(See  reverse  side) 


79  09 


0 27 


DD  , 1473  COITION  OF  I NOV  ••  IS  OBSOLETE 

S/N  0102- LF-01 4-6601 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  <•*>>•"  »olo  BntotoO) 


00 $ $00 


SECURITY  CLASSIFICATION  of  THIS  PAGE  CHTiwi  Data  Entered) 


Aeronautical  Systems  Division 
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One  of  the  problems  involved  in  computer  programs  for 
vulnerability  assessment  is  that  of  rapidly  sorting  and  arrang- 
ing large  sets  of  data.  Two  sorting  algorithms,  designated 
PRIME  and  PDQ,  have  been  developed  at  ASD  to  more  effi- 
ciently perform  this  function  in  vulnerability  programs  such 
as  SESTEM  and  FASTGEN  II.  The  results  are  compared  to 
those  obtained  with  three  other  algorithms,  SHELLSORT, 
TREESORT3,  and  SINGLETON.  The  newly  developed  sorts 
are  shown  to  be  significantly  faster  on  the  ASD  CDC  6600 
computer  than  the  existing  sorts.  When  used  in  an  ASD  mis- 
sile endgame  model  SESTEM,  the  average  run  time  was 
reduced  by  20  to  25%.  Program  listings,  flow  charts,  and 
typical  output  data  are  presented.  ^ 
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INTRODUCTION 
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One  of  the  steps  in  most  automated  vulnerability  assessments  involves  the  rapid  sorting 
of  large  strings  of  data.  In  the  ASD  (Aeronautical  Systems  Division)  missile  endgame  model 
SESTEM,1  for  example,  the  aircraft  components  struck  by  the  expanding  fragment  spray- 
band  are  calculated  and  stored,  and  then  must  be  sorted.  The  sorted  data  are  then  used  to 
compute  the  components  struck  in  order  of  time  intercepted  and  the  resulting  aircraft  prob- 
ability of  kill.  Similar  data  string  storage  and  sort  problems  are  involved  in  other  vulnerabil- 
ity assessment  computer  programs  such  as  FASTGEN  II.2  In  SESTEM,  about  40 % of  the 
computation  time  is  typically  involved  with  the  sorting  and  fragment  vulnerability 
computations. 

For  these  reasons,  the  development  of  efficient  algorithms  to  sort  these  strings  of  data 
are  of  considerable  importance.  Since  large  portions  of  the  total  computer  run  times  are  typ- 
ically involved  in  performing  this  sorting,  any  improvements  in  them  will  show  direct  pay- 
offs in  terms  of  decreased  run  times,  more  rapid  turnaround,  greater  program  efficiencies, 
and  decreased  costs.  The  sorting  algorithms  in  use  in  the  vulnerability  assessment  programs 
were  evaluated,  and  two  new  and  more  efficient  programs  were  written.  When  the  sorting 
algorithm  was  inserted  in  a new  version  of  the  SESTEM  program,  the  average  time  for  a run 
was  decreased  by  20  to  25%.  These  programs,  and  comparisons  with  some  existing  programs 
are  presented  in  this  report. 


PROGRAM  DEVELOPMENT 


BACKGROUND 

The  availability  of  efficient  general  purpose  sort  algorithms  and  the  sharp  reduction  of 
cost  per  computation  of  present  generation  computers  has  reduced  interest  in  research  into 
sort  algorithms  from  earlier  levels.  Aside  from  the  intrinsic  challenge  presented  by  the 
sorting  problem,  optimization  of  existing  and  developing  applications  programs  written  in 
FORTRAN  IV,  or  similar  level  languages,  indicate  a still  existing  need  for  compact, 
efficient,  in-line  sort  algorithms  that  are  readily  adaptable  to  specialized  ends. 

One  approach  to  meeting  this  need  is  to  devise  a partial  sort,  an  efficient  algorithm  by 
which  a string  of  numbers  is  nearly  sorted;  and  complete  the  sort  by  a method  such  as 
binary  insertion,  which  can  take  advantage  of  a high  degree  of  order  in  a number  string. 

* Aeronautical  Systems  Division.  “SESTEM  Missile  Endgame  Model",  by  G.B.  Bennett,  Wright-Patterson  Al  B,  OH, 
ASD.  September  1977.  (ASD/XRH  memo.) 

2 Aeronautical  Systems  Division.  FASTGEN  It  Target  Description  Computer  Program,  by  D.  Cudney. 
Wright-Patterson  AEB.OH,  ASD.  March  1978.  (Report  ASD-TR-77-24.) 
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PARTIAL  SORTING 

The  basic  strategy  in  producing  partially  sorted  strings  is  to  successively  partition  a set 
of  numbers  by  exchange  comparisons.  Let  NUM  be  a string  of  numbers  of  length  NO  equal 
to  an  integral  power  of  two.  The  set  of  NO  numbers  is  divided  into  two  subsets  of  equal  size, 
and  the  elements  of  one  subset  are  compared  (and  exchanged  when  necessary)  to  the  ele- 
ments of  the  other  subset  in  such  a way  that  every  element  is  involved  in  a comparison  once 
only.  As  a worst  case,  after  any  comparison  of  sets  as  described,  25%  of  the  elements  less 
than  or  equal  to  the  median  value  will  be  located  below  the  median.  Successive  set  divisions 
and  comparisons  continue  the  process,  removing  extreme  values  from  the  middle  of  the 
array  and  enabling  further  distribution  to  take  place,  until  there  are  NO  subsets  of  1 ele- 
ment each.  If  the  NO  elements  are  distinct  and  represented  by  1,2,...,  NO,  Figure  1 
represents  the  process  described. 


PDQ  ALGORITHM 

The  successive  division  of  subsets  into  equal  subsets,  as  illustrated,  will  cause  the 
1 element  to  migrate  to  the  proper  postion  in  the  array ; but  in  the  worst  case,  the  2 element 
will  migrate  toward  the  middle  position  of  the  array.  The  2 element  can  be  forced  left  by 
dividing  the  subsets  unequally,  as  shown  in  Figure  2.  Let  the  kth  interval  of  comparison  be 
defined  as  the  integer  part  of  1^  = NO*F^,  where  0.<  F<  1 . A code  to  achieve  the  partition- 
ing is  given  in  Listing  1.  F = .8  seems  to  produce  a complete  sort  on  a uniform  random 
distribution;  however,  the  author  has  not  been  able  to  characterize  the  optimal  value  for  F 
for  a particular  given  distribution.  Observation  of  empirical  test  results  suggests  that  certain 
structured  distributions  are  difficult  to  sort;  therefore,  the  PDQ  algorithm  may  be  useful  as 
a test  of  randomness. 


SORT  CODE  TESTING  ALGORITHM 


DISTRIBUTIONS 

The  sort  code  was  tested  by  a control  program  (Listings  2A  to  2D)  that  generates 
strings  of  a specified  length  from  six  different  distributions; 

1 . Sorted  arrays:  the  elements  of  the  sequence  1,2,3,...,  NO. 

2.  Arrays  in  reverse  order:  the  elements  of  the  sequence  NO,  NO-1  ,...,2,1. 

3.  Random  arrays:  numbers  from  a uniform  distribution  over  the  interval  (0.  , 1.). 
The  random  numbers  are  multiplied  by  100,000  to  minimize  duplicate  numbers  when  con- 
verted to  integers. 
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4.  Arrays  almost  in  sort:  the  sequence,  1,2,...,  NO  with  a specified  number  of  ele- 
ments, chosen  at  random,  set  to  values  taken  from  a uniform  random  distribution. 

5.  Arrays  of  equal  length  sorted  blocks:  this  is  the  distribution  described  in  (3)  sorted 
in  successive  segments  of  a specified  length. 

6.  Constant  value  arrays:  a sequence  of  numbers  of  equal  value.  These  distributions 
were  selected  with  a view  to  test  algorithm  behavior  on  distributions  that  might  be  encoun- 
tered in  practical  applications,  such  as  (3)  and  (4),  or  to  demonstrate  unusual  characteristics. 

Elapsed  time  to  sort  is  measured,  and  a count  of  departures  from  a monotonic  sequence 
is  made.  If  this  count  is  greater  than  0,  an  error  message  is  printed. 


TESTING  ENVIRONMENT 

All  test  runs  were  conducted  on  a CDC  CYBER  74,  with  the  same  level  of  optimization 
(OPT  = 2).  Observations  of  repeated  runs  on  the  CYBER  74,  operating  in  a time-sharing 
mode,  suggest  that  time  measurements  can  vary  about  20%;  howeverrby  specifying  large 
arrays  to  sort  (90K),  the  job  will  be  made  to  run  on  the  machine  in  a more  dedicated  con- 
figuration. In  this  dedicated  mode,  elapsed  times  are  highly  reproducible.  Sort  times  are 
sensitive  to  the  values  taken  for  F,  as  shown  in  Figure  3. 

Machine  Dependancy 


An  apparent  anomally  is  that  it  takes  more  time  to  sort  a sorted  string  than  to  sort  a 
string  that  has  an  initial  uniform  random  distribution.  This  seems  to  indicate  that  it  takes 
more  time  to  execute  a branch  instruction  than  the  three  arithmetic  replacement  statements 
involved  in  the  number  exchange.  This  is  true  in  the  aggregate  for  the  repetitive  execution  of 
the  code  in  the  DO  loop  used  to  compare  and  exchange  the  elements  of  the  NUM  array. 

The  CYBER  74  is  a stack  machine.  Up  to  seven  words  of  packed  instructions  contain- 
ing up  to  28  instructions  can  be  retained  in  registers  constituting  an  instruction  stack.  This 
device  can  increase  instruction  execution  speed  by  reducing  memory  references;  however,  a 
forward  branch  in  the  stack  “voids  the  stack”,  therefore  is  an  expensive  operation.  For  this 
reason,  PDQ  will  sort  faster  if  the  “less  than  equal”  test  is  replaced  by  a “less  than”  test  in 
the  exchange  algorithm.  This  increased  speed  is  demonstrated  in  the  decreased  times  required 
to  sort  a constant  value  array,  as  compared  to  the  time  required  to  sort  an  array  already  in 
sort. 


Markedly  different  results  may  be  expected  on  some  other  computer  systems.  The 
expected  time  relationship  may  be  obtained  by  replacing  the  DO  loop  with  an  IF  loop, 
where  the  branch  will  be  in  a backward  direction  in  the  instruction  stack  (Listing  3).  Unfor- 
tunately, code  optimization  is  not  as  intensive  now,  and  overall  sort  times  are  increased. 
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Figure  1.  Successive  Partition  of  Sets  into  Equal  Subsets. 
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Figure  2.  PDQ  Partitioning  Process. 
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LISTING  1 
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TIME  IN  SECONDS  AND  FPROR  COUNTS  FOR  TWO 
SUCCESSIVE  POO  SORTS  CN  A UNIFORM  RANDOM 
DISTRIBUTION  OF  HOBO  NUMBERS 
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Figure  3.  Elapeed  Time  as  a Function  of  F for  the  PDQ  Partial 
Sort  (Listing  1). 
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The  simplicity  of  the  PDQ  code  lends  it  to  numerous  variations;  for  instance,  PDQ 
modules  can  be  stacked  with  different  values  for  F in  each  module.  Changing  the  direction 
of  the  comparisons  on  alternate  passes  causes  sorting  to  occur  in  fewer  passes  but  at  the  cost 
of  increased  complexity  (Listing  4).  Listing  5 is  a format  of  PDQ  using  a partitioning  scheme 
involving  the  logarithm  to  the  base  12. 


C 

10 


15 


2C 

C 


gMTO,?W!No 


\ 


IF  1 1 

< 


= NO 
x A*F 

* A 
.LE.  0) 

* NO  - I 


RETURN 


DO  15  J * l.K 

XFINUMtut.LE.NUMtXtjn  GO  TO  15 


MAX 


* NUM  t J) 


NUM  ( J)  = N(JM(I*  J) 
KUMjI*J>=  MAX 
CONTINUE 


If 1 1 .le. 
< * 

L s 

CO  20  J * 


A*F 

A 

0) 

NO  - X 
K ♦ 1 

l'K- 


RETURN 


IF(NUM(L) .LC.NUMt I 
MAX  « NUM ( LI 
NUMt  LI  - NUM(I»LI 


*LD  GO  TO  20 


NUM  tl*L>*  MAX 
CONTINUE 
GO  TO  ll 

END 


POQ  ALTERNATE 


LISTING  k 


8 


JTCG/  AS-78-V-004 


10 


15 


ay  ««f 


POO  LCC12 


% 


P 

A 

irti  « 

K 

00  15 


LE. 
J = 


NUH2. 

IMNUN1  , 

NUN ( J)  i 
NUMtI»J>=  MINI 
CONTINUE 
GO  TO  10 

END 


NO 

ALOG16(AI/l. 07918 

P - • 1 

12..-P 

5 I 

NO  - I 

l.K 

NUN ( J » 

NJNII» J> 

LT.  NUN 2 ) 

NUH2 


RETURN 
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LISTING  5 


The  original  intention  of  the  sort  strategy  is  achieved,  with  the  assurance  of  a sort  by 
following  PDQ  with  a sort  by  direct  insertion.  The  insertion  code  is  derived  from  PDQ  by 
the  addition  of  a few  lines  of  code  (Listing  6). 

Insertion  means  the  addition  of  numbers  to  an  existing  sorted  string  (which  initially 
may  be  of  length  1 ) by  inserting  the  new  element  into,  or  at  the  ends  of  the  sorted  string  to 
form  a new  sorted  string  of  increased  length.  For  any  distribution,  the  sort  is  completed  in 
one  pass.  An  INSERTION  sort  is  to  be  distinguished  from  a BUBBLE  sort  where  successive 
elements  are  selected  and  added  to  one  end  of  a string,  initially  of  length  0.  The  BUBBLE 
sort  may  require  up  to  NO-1  passes  to  complete  the  sort.  Tests  show  the  INSERTION  sort 
to  be  more  efficient  than  the  BUBBLE  sort  (Listing  7). 

If  the  features  of  the  PDQ  sort  and  the  INSERTION  sort  are  combined,  a particularly 
efficient  algorithm  resembling  the  SHELL  sort  is  obtained  (Listing  8).  The  specification  of 
an  optimal  sequence  of  intervals  to  control  partitioning  is  a difficult  task;  however,  if  a geo- 
metric sequence  is  assumed,  an  F can  empirically  be  found  which  will  yield  minimum  sort 
times  for  a given  distribution  (Figure  4).  If  two  or  more  values  of  F produce  minimal  elapsed 
times,  the  smallest  of  these  values  should  be  selected  to  reduce  sort  times  on  sorted  or 
nearly  sorted  strings. 

The  SHELL  sort  seems  to  have  been  intended  as  a type  of  merge  algorithm.3  The  term 
“merge”  may  be  used  in  several  ways: 

1 . The  combination  of  two  or  more  sorted  strings  into  a resultant  sorted  string  irre- 
spective of  any  particular  algorithm. 

2.  A specific  method  by  which  sorted  strings  can  be  combined  efficiently  into  result- 
ant sorted  strings. 


3Gotlieb,  C.C.  Sorting  on  Computers  Communications.  ACM  6 (May  1963),  p.  194. 
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Figure  4.  Elapsed  Time  as  a Function  of  F for  Distribution  1 (Listing  8). 
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Generally,  the  SHELL  sort  does  not  qualify  as  a merge  algorithm,  but  it  can  be  made 
to  operate  as  such  under  definition  (1 ) by  selecting  F = .5  and  the  first  interval  as  the  largest 
power  of  2 less  than  NO. 

This  choice  of  F turns  out  to  be  one  of  the  worst  possible.  To  see  why,  it  may  be  con- 
structive to  consider  this  procedure  as  a type  of  distribution  sort.  Consider  the  sequence  of 
16  numbers,  1,  2,  . . 16  where  the  integers  represent  the  position  in  an  array  of  numbers 

to  be  sorted.  Take  the  first  interval  of  comparison,  I|  equal  to  8 and  F = .5,  then  the  follow- 
ing illustration  can  be  used: 

1 2 3 4 5 6 7 8 Pass  I 

9 10  11  12  13  14  15  16  1 1 = 8 

to  suggest  that  the  array  to  be  sorted  has  been  divided  into  eight  subsequences  represented 
by  the  columns.  The  numbers  in  each  column  are  to  be  sorted  in  ascending  order  from  top 
to  bottom;  i.e.,  the  number  in  position  1 is  to  be  less  than  or  equal  to  the  number  in  posi- 
tion 9,  etc. 


The  following  illustration  represents  the  next  pass: 

1 

2 

3 

4 

Pass  2 

5 

6 

7 

8 

I2  = 4 

9 

10 

11 

12 

13 

14 

15 

16 

Note  that:  (1 ) the  subsequences  are  reduced  in  number  in  proportion  to  F,  (2)  the  length  of 
the  subsequences  is  increased  in  inverse  proportion  to  F,  and  (3)  the  elements  of  a given 
column  are  composed  alternately  of  the  n and  n + lj/2  columns  of  the  previous  pass, 
n = 1,  2,  3,  Ij/2.  Consequently,  there  is  a high  degree  of  order  in  any  column  and  an  ele- 
ment is  likely  to  be  close  to  its  sorted  position  in  a subsequence,  and  an  INSERTION  sort  or 
BUBBLE  sort  applied  to  a column  may  be  expected  to  operate  faster  than  on  a random 
distribution  of  the  same  length.  The  worst  case  in  this  example  occurs  when  the  even  and 
odd  columns  represent  disjoint  ranges,  and  the  set  of  numbers  from  the  even  columns 
contains  the  smaller  numbers.  Insofar  as  the  median  of  each  column  represents  the  median 
of  the  entire  distribution,  subsequent  passes  may  be  expected  to  produce  increasingly 
well-ordered  subsequences  representative  of  the  entire  sequence;  therefore,  needing  a 
minimum  number  of  comparisons  and  exchanges  to  sort.  However,  certain  difficulties  can 
arise. 


In  the  following  example,  the  array  notation  of  the  previous  illustration  will  be  re- 
tained, but  now  the  integers  represent  the  actual  numbers  in  the  array  that  are  to  be  sorted. 
Given  the  sequence  3,  5,  6,  7,  2,  8,9,  10,4,  11,  12,  13,  1,  14,  15,  16,  let  F = .5  and  Ij  = 8, 
we  have : 


3 5 6 7 


8 9 10  Pass  1 


4 11  12  13 


14  15  16  1 1 = 8 


r« 


- * ■ 


2 
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AH  columns  but  the  fifth  are  in  sort;  therefore,  the  1 and  2 elements  are  exchanged.  This 
could  be  called  an  unfavorable  exchange,  as  it  improves  the  order  of  the  array  very  little; 
both  the  1 and  the  2 should  be  in  the  first  row.  After  the  exchange,  pass  2 can  be 
represented : 


3 

5 

6 

7 

Pass  2 

1 

8 

9 

10 

I2  = 4 

4 

11 

12 

13 

2 

14 

15 

16 

The  columns  are  well-ordered  except 

column  1 

. After  the  sort  we  have: 

1 

5 

6 

7 

2 

8 

9 

10 

3 

11 

12 

13 

4 

14 

15 

16 

Column  1 has  a median  atypical  of  the  array,  and  most  of  its  elements  are  far  from  their 
final  position.  It  will  take  a large  number  of  small  steps  to  move  them  into  postion  in  later 
passes.  The  elements  in  column  1 after  row  1 may  be  thought  of  as  blocking  elements  be- 
cause they  inhibit  efficient  distribution  of  the  array  elements  in  the  early  stages  of  the  sort 
when  an  element  can  proceed  toward  its  destination  by  large  steps.  To  follow  the  example 
further; 

1 

5 

Pass  3 1 

, 5 

6 

7 

I3  = 2 2 

7 

2 

8 

3 

8 

9 

10 

4 

10 

3 

11 

6 

11 

12 

13 

9 

13 

4 

14 

12 

14 

15 

16 

15 

16 

The  last  pass  with  I4  = 1 will  require  a large  number  of  exchanges  to  complete  the  sort. 
The  collection  of  blocking  elements  can  be  prevented  by  a number  of  methods,  such  as  sort- 
ing up  a diagonal  from  left  to  right,  following  a column  sort.  But  most  of  the  comparisons 
will  not  result  in  an  exchange  due  to  a high  degree  of  order  produced  by  the  previous  pass, 
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and  those  exchanges  that  do  result  are  bought  at  the  expense  of  an  extra  pass.  Similar  objec- 
tions may  be  raised  to  other  modifications  to  the  algorithm  to  eliminate  the  blocking  ele- 
ments. Examination  of  the  problem  reveals  that  the  accumulation  of  blocking  elements 
from  pass  to  pass  is  particularly  severe  if  the  rows  in  the  array  are  divided  into  an  integral 
number  of  rows  in  the  next  pass;  therefore,  .25  and  .5  are  bad  values  for  F.  This  can  be 
verified  readily  by  rewriting  the  sort  algorithm  to  include  counts  of  primary  exchanges  and 
secondary  exchanges  on  a pass-by-pass  basis,  and  printing  these  values  (as  well  as  their  cumu- 
lative totals),  for  selected  values  of  F (Figures  5 and  6).  It  is  interesting  to  note  that  the  mini- 
mum sort  times  do  not  correspond  to  the  values  of  F that  result  in  minimum  exchanges, 
because  as  F becomes  larger,  there  is  an  increase  in  overhead  associated  with  the  increase  in 
the  number  of  passes.  The  elapsed  times  shown  are  distorted  by  the  code  modification 
required  for  the  accumulation  and  printing  of  the  statistics. 

A study  of  the  exchange  counts  suggests  that  the  blocking  phenomena  come  into  play 
whenever  there  is  a series  of  intervals  that  have  factors  in  common,  and  that  the  problem 
becomes  more  severe  as  F decreases;  therefore,  the  intervals  should  be  relatively  prime.  Con- 
sider the  first  column  of  the  array  on  the  mth  pass.  The  location  of  a blocking  element  in 
the  first  column  can  be  expressed  as  1 + M*Pm.  On  the  next  pass,  if  that  same  element  is  to 
appear  in  the  same  column,  it  must  have  a location  expressible  as  I + N*Pn,  where  Pn  < Pm. 
Then  N*Pn  = M*Pm,  and  N = Pm,  M = Pn  if  Pm.  Pn  are  relatively  prime.  But  these  values  for 
M and  N are  impossible  during  the  earlier  passes  of  the  sort  for  large  arrays  and  the  range  of 
F considered  here.  To  test  this  hypothesis,  sequences  of  relatively  prime  numbers  corre- 
sponding to  the  geometric  sequences  associated  to  values  for  F (F  = .2  to  F = .46)  were  used 
in  the  sort  algorithm  (Listing  9)  to  conduct  elapsed  time  tests  (Figure  7).  The  tests  seem  to 
confirm  the  hypothesis;  however,  the  use  of  prime  sequences  is  only  a partial  solution.  A 
blocking  element  may  be  replaced  in  a column  with  another  element  that  also  serves  as  a 
blocking  element.  This  situation  is  likely  to  occur  when  sorting  arrays  of  equal  length  sorted 
blocks.  Counts  of  exchanges  and  comparisons  for  the  prime  sequences  are  given  in  Figure  8. 


PRIME  SORTING  ALGORITHM 


The  prime  sequence  corresponding  to  F = .3  was  selected  with  a view  to  minimum  com- 
bined elapsed  times  for  both  random  and  sorted  distributions  for  the  algorithm  PRIME 
SORT  (Listing  10).  This  algorithm  also  incorporates  a change  that  significantly  improves  the 
efficiency  of  exchanges,  especially  secondary  exchanges. 

The  improvement  in  execution  speed  was  largely  lost  when  the  data  was  passed  through 
the  CALL  list.  This  illustrates  a compilation  problem  that  affects  the  various  algorithms 
given  here  to  differing  extents;  the  PDQ  codes  were  degraded  least  by  passing  data  through 
the  CALL  list.  Efficient  compilation  is  important  because  the  PDQ,  PRIME,  and  DISTRI- 
BUTION sort  codes  given  here  derive  their  speed  from  compact  code  that  requires  a mini- 
mum of  instructions  to  execute,  and  from  improved  partitioning  schemes.  Full  realization  of 
the  potential  of  the  code  requires  effective  register  assignment  to  indices,  etc.,  but  variables 
passed  through  CALL  lists  inhibit  optimal  compilation.  Generally,  an  improvement  in  per- 
formance may  be  expected  when  data  is  passed  to  a subroutine  through  COMMON  storage 
(Figure  9). 
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Figure  6.  Cumulative  Count*  of  Exchanges  and  Comparisons  as  a 
Function  of  F for  the  Distribution  1 1 Algorithm. 
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Figure  7.  Elapsed  Time  for  Prime  Sequences,  F * .2  to 
F * .46. 
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Figure  8.  Cumulative  Counts  of  Exchanges  and  Comparisons  for 
Prime  Sequences,  F * .2  to  F * .46. 
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Assembly  language  coding  can  minimize  these  problems,  and  candidates  for  that  pur- 
pose are  given  in  Listings  1 1 and  12.  An  algorithm  suitable  for  in-line  code  is  given  in  List- 
ing 13.  The  effect  of  a prime  sequence  is  approximated  by  making  all  intervals  odd,  etc. 
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Distribution  2 (Listing  13). 
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For  each  primary  comparison,  there  is  a chance  of  a primary  exchange  which  will  be 
followed  generally  by  a secondary  comparison,  etc.  Figure  6 indicates  that  secondary  ex- 
changes will  be  substantial  for  useful  values  of  F;  therefore,  the  expression  for  the  compari- 
sons of  the  DISTRIBUTION  or  PRIME  sort  will  be  that  for  PDQ  sort,  plus  a series  of  terms 
involving  probabilities  which  represent  the  various  orders  of  exchanges  (primary,  secondary, 
etc.).  Sort  times  for  large  NO  may  be  expected  to  be  proportional  to  Log  NO  for  the  PDQ 
sort,  and  to  increase  faster  than  a Log  NO  rate  for  PRIME  sort. 


PERFORMANCE  EVALUATION 


For  a comparative  evaluation  of  sort  performance,  some  establised  sort  algorithms 
published  in  the  COMMUNICATIONS  of  the  ACM4'S * 7  were  used.  These  algorithms  were 
adapted  for  the  sake  of  uniform  notation  and  style.  In  addition,  the  subroutine  SIFTUPof 
TREESORT3  was  coded  in-line  to  reduce  the  substantial  overhead  involved  in  the  frequent 
calls  to  this  procedure.  The  sort  designated  SINGLETON  is  one  of  the  faster,  more  stable 
members  of  the  QUICKSORT  family. 


4 

Singleton,  Richard  C.  An  efficient  algorithm  for  sorting  with  minimal  storage.  Communications  ACM 
12  (March  1969),  p.  185. 

SLoeser,  Rudolf.  Some  performance  tests  of  “QUICKSORT”  and  decendants.  Communications  ACM 
17  (March  1974),  p.  143. 

®Boothroyd,  J.  Algorithm  201,  SHELLSORT.  Communications  ACM  6,  8 (August  1963),  p.  445. 

7Eloyd,  R.  W.  Algorithm  245,  TREESORT3.  Communications  ACM  7,  12  (December  1964),  p.  701. 
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An  attempt  to  evaluate  the  efficiency  of  a procedure  by  frequency  counts  of  critical 
parameters  is  not  entirely  satisfactory,  even  for  differing  versions  of  that  procedure  on  the 
same  machine.  The  particular  form  of  a procedure  is  very  important.  Machine  independent 
comparisons  between  differing  procedures  are  even  more  difficult.  Some  of  the  results  were 
omitted  from  the  graphical  sort  performance  data  when  interference  between  curves  obscured 
the  comparisons.  Each  curve  has  a label  that  refers  to  the  listing  of  the  code  used  in  generat- 
ing the  data,  (e.g.,  PDQ  I L18  means  PDQ  INSERT  #2,  Listing  18.)  The  test  results  indicate 
that  the  PRIME  and  DISTRIBUTION  sorts  compare  favorably  overall  to  the  QUICKSORT ; 
and  in  the  case  of  random  distributions,  the  results  increasingly  favor  the  PDQ,  PRIME,  and 
DISTRIBUTION  sorts  as  the  array  size  decreases.  The  DISTRIBUTION  sort  is  a compact 
and  efficient  sort  suitable  for  in-line  code  applications  because  it  is  generally  understandable; 
therefore,  it  may  be  modified  easily  for  particular  uses. 
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SUBROUTINE  SORT (NUHtNC) 

(ACM  ALGORITHM  3A7)  BY  RICHARD  C.  SINGLETON 
. COMMUNICATIONS  OF  THE  ACM  - VOLUME  12  / NUMBER  3 / MARCH.  1969 
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Aeronautical  Systems  Division  (AFSC) 

Wright-Pat ter son  AFB,  OH  45433 
Attn:  ASD/ENESH  (P.  T.  Marth) 

Attn:  ASD/ENFTV  (D.  J.  Wallick)  (2  copies) 
Attn:  ASD/XROL  (F.  Campanile) 

Attn:  ASD/XROL  (R.  K.  Frick) 

Attn:  ASD/XROT  (G.  B.  Bennett) 

Attn:  ASD/XRE  (S.  E.  Tate) 

Attn:  ASD/YPES  (C.  Gebhard) 

Aerospace  Medical  Research  Laboratories 

AMRL/MEA 
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Wright-Patterson  AFB,  OH  45433 

Attn:  AMRL/MEA  (CAPT  G.  J.  Valentino) 

Aerospace  Rescue  and  Recovery  Service  (MAC) 
Scott  AFB,  IL  62225 

Attn:  ARRS/DOQ  (CAPT  J.  J.  Draham,  Jr.) 

Air  Force  Acquisition  Logistics  Division 
Wright-Patterson  AFB,  OH  45433 

Attn:  AFALD/PTEA  (MAJ  D.  Waltman) 

Air  Force  Avionics  Laboratory 
Wright-Patterson  AFB,  OH  45433 
Attn:  AFAL/WRA-1  (E.  Leaphart) 

Attn:  AFAL/WRP  (W.  F.  Bahret) 

Air  Force  Flight  Dynamics  Laboratory 
Wright-Patterson  AFB,  OH  45433 

Attn:  AFFDL/FES  (CDIC)  (2  copies) 

Attn:  AFFDL/FES  (C.  W.  Harris) 

Attn:  AFFDL/FES  (J.  Hodges) 

Attn:  AFFDL/FES  (R.  W.  Lauzze) 

Attn:  AFFDL/FES  (D.  W.  Voyls) 

Attn:  AFFDL/FGL 

Attn:  AFFDL/TST  (Library) 

Air  Force  Logistic  Command 
Wright-Patterson  AFB,  OH  45433 
Attn:  AFLC/LOE  (Commander) 

Air  Force  Systems  Command 
Andrews  AFB,  DC  20334 

Attn:  AFSC/DLCAA  (P.  L.  Sandler) 
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Air  Force  Test  and  Evaluation  Center 
Kirtland  AFB,  NM  87115 

Attn:  AFTEC/OAR  (K.  Campbell) 

Attn:  AFTEC/OAR  (R.  Blankert) 

Air  Force  Weapons  Laboratory 
Kirtland  AFB,  NM  87117 

Attn:  AFWL/PGV  (CAPT  J.  K.  Carson) 

Attn:  AFTEC/OA  (MAJ  H.  Rede) 

Applied  Technology  Laboratory 

Army  Research  & Technology  Laboratory  (AVRADCOM) 
Ft.  Eustis,  VA  23604 

Attn:  DAVDL-ATL-AL  (Mr.  Merritt) 

Attn:  DAVDL-ATL-ASV  (S.  Pociluyko) 

Attn:  DAVDL-ATL-ASV  (H.  W.  Holland) 

Attn:  DAVDL-ATL-ASV  (J.  T.  Robinson) 

Armament  Development  and  Test  Center 
Eglin  AFB,  FL  32542 

Attn:  ADTC/DLODL  (Technical  Library) 

Attn:  ADTC/XR  (C.  T.  Maney) 

Attn:  ADTC/XRSP  (E.  Blair) 

Armament  Research  and  Development  Command 
Dover  Base 
Dover,  NJ  07801 

Attn:  DRDAR-LCS  (S.  K.  Einbinder) 

Army  Aviation  Research  & Development  Command 

P.0.  Box  209 

St.  Louis,  MO  63166 

Attn:  DRCPM-ASE-TM  (MAJ  Schwend)  (2  copies) 

Army  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground,  MD  21005 
Attn:  DRDAR-BLV  (D.  W.  Mowrer) 

Army  Foreign  Science  and  Technology  Center 
220  Seventh  St.,  NE 
Charlottesville,  VA  22901 

Attn:  DRXST-BA3  (E.  R.  Mclnturff) 

Attn:  DRXST-CA2  (J.  M Blake)  (2  copies) 

Army  Materials  and  Mechanics  Research  Center 
Watertown,  MA  02172 

Attn:  DRXMR-ER  (F.  C.  Quigley) 

Attn:  DRXMR-EM  (C.  F.  Hickey,  Jr.) 

Attn:  DRXMR-PL  (M.  M.  Murphy)  (2  copies) 
Attn:  DRXMR-XC  (E.  S.  Wright) 
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Army  Materiel  Systems  Analysis  Activity 
Aberdeen  Proving  Ground,  MD  21005 
Attn:  DRXSY-AA  (Director) 

Attn:  DRXSY-AD  (H.  X.  Peaker) 

Attn:  DRXSY-J  (J.  J.  McCaitrhy) 

Attn:  DRXSY-S  (J.  R.  LindenautM 

Chief  of  Naval  Operations,  Office  of  the 
Room  5C735,  Pentagon 
Washington,  DC  20350 

Attn:  OPNAV-982E3L  (LT  COL  W.  A.  Allanson) 

David  W.  Taylor  Naval  Ship  R&D  Center 
Carderock  Laboratory 
Bethesda,  MD  20084 

Attn:  Code  1740.2  (F.  J.  Fisch) 

Attn:  Code  1740.2  (F.  Hackett) 

Attn:  Code  1740.4  (M.  L.  Salive) 

Attn:  Code  522  (Commander)  (4  copies) 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Blvd. 

Arlington,  VA  22209 

Attn:  S.  Zakanycz  DARPA/STO 

Defense  Documentation  Center 
Cameron  Station,  Bldg.  5 
Alexandria,  VA  22314 

Attn:  DDC-TCA  (2  copies) 

Defense  Systems  Management  College 
Ft.  Belvoir,  VA  22060 
Attn:  W.  Schmidt 

Department  of  Transportation  - FAA 
2100  Second  St.,  SW,  Rm  1400C 
Washington,  DC  20591 

Attn:  ARD-520  (R.  A.  Kirsch) 

Deputy  Chief  of  Staff  (AIR) 

Marine  Corps  Headquarters 
Washington,  DC  20380 

Attn:  APW-22  (MAJ  K.  Van  Esselstyn) 

ERADCOM 

Fort  Monmouth,  NJ  07703 

Attn:  DAVAA-I  (S.  Zywotow,  Avionics  R&D  Activity) 

Attn:  DELEW-P  (R.  F.  Giordano,  Electronic  Warfare  Laboratory) 
Attn:  DELSD-E  (C.  Goldy) 

Attn:  DRSEL-GG-TD  (Commander) 
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Foreign  Technology  Division  (AFSC) 

Wright-Pat terson  AFB,  OH  45433 
Attn:  FTD/NICD  (2  copies) 

HQ  Air  Logistics  Command 
McClellan  AFB,  CA  95652 

Attn:  SM/MMSRBC  (D.  E.  Snider) 

Joint  Strategic  Target  Planning  Staff 
Offutt  AFB,  NB  68113 

Attn:  JSTPS/JPTB  (MAJ  C.  0.  Cox) 

Attn:  NRI/STINFO  (Library) 

NASA  - Ames  Research  Center 

Army  Research  and  Technology  Laboratories,  Headquarters 

Mail  Stop  207-5 

Moffett  Field,  CA  94035 

Attn:  DAVDL-AS  (V.  L.  J.  Di  Rito) 

NASA  - Johnson  Space  Center 
Houston,  TX  77058 

Attn:  EC  (F.  S.  Dawn) 

Naval  Air  Development  Center 
Warminster,  PA  18974 

Attn:  Code  097  (MAJ  W.  Boeck) 

Attn:  Code  701  (B.  Vafkos) 

Attn:  Code  2012  (C.  E.  Murrow) 

Attn:  Code  2012  (M.  C.  Mitchell) 

Attn:  Code  2012  (R.  H.  Beliveau) 

Attn:  Code  3023  (L.  M.  Rakszavski) 

Attn:  Code  6083  (S.  L.  Huang) 

Attn:  Code  601A  (R.  A.  Ritter) 

Naval  Air  Propulsion  Center 
P.0.  Box  7176 
Trenton,  NJ  08628 

Attn:  PE3  (D.  Wysocki) 

Naval  Air  Systems 
Airtevron  One 
Patuxent  River,  MD  20653 
Attn:  LT  R.  N.  Freedman 
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Naval  Air  Systems  Command 
Washington,  DC  20361 

Attn:  AIR-03PA4  (T.  S.  Momiyama) 

Attn:  A1R-330B  (E.  A.  Lichtman) 

Attn:  AIR- 350  (E.  M.  Fisher) 

Attn:  AIR-503W1  (E.  A.  Thibault) 

Attn:  AIR-52014  (L.  Sztan) 

Attn:  AIR-5184  (D.  Atkinson)  (2  copies) 
Attn:  AIR- 5184 J (MAJ  R.  A.  Horton) 

Attn:  AIR-5303 

Attn:  AIR-530313  (R.  D.  Hume) 

Attn:  AIR-53051A  (P.  Kicos) 

Attn:  AIR- 532 3A 

Attn:  AIR-5323K  (K.  Chang) 

Attn:  AIR-5323Z  (S.  Englander) 

Attn:  AIR-620B1  (LCDR  K.  K.  Miles) 

Attn:  AIR- 954  (Tech.  Library) 

Attn:  PMA-2692A1  (R.  W.  Wills) 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Attn:  Code  67BP  (R.  E.  Ball) 

Attn:  Library 

Naval  Research  Laboratory 
4555  Overlook  Ave. , SW 
Washington,  DC  20375 

Attn:  Code  1409  (J.  M.  MacCallum) 

Attn:  Code  5730  (E.  E.  Koos) 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

Attn:  SEA-6543  (F.  W.  Sieve) 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22448 

Attn:  CK-2301  (J.  E.  Mitchell) 

Attn:  CN-61  (J.  S.  Nerrie) 

Attn:  G-10  (J.  E.  Ball) 

Attn:  G-10  (S.  Hock) 

Attn:  G-10  (F.  J.  Petranka) 

Attn:  Library 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 
Attn:  CN-13 
Attn:  WX-21  (Library) 
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Naval  War  College 
Newport,  RI  02840 

Attn:  Code  E-lll  (President) 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Attn:  Code  317  (M.  H.  Keith) 

Attn:  Code  3181  (C.  Padgett)  (2  copies) 

Attn:  Code  3183  (C.  Driussi) 

Naval  Weapons  Engineering  Support  Activity 
Systems  Analysis  Dept. 

Bldg.  210-2  (ESA-19) 

Washington  Navy  Yard 
Washington,  DC  20374 

Attn:  Code  ESA-1923  (C.  W.  Stokes  III)  (2  copies) 

Attn:  Code  11621  (J.  Stasko) 

Naval  Weapons  Support  Center 
Crane,  IN  47522 

Attn:  Code  502  (N.  L.  Papke) 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217 

Attn:  0NR  474  (N.  Perrone) 

Pacific  Missile  Test  Center 
Point  Mugu,  CA  93042 

Attn:  Code  4253-3  (Naval  Air  Station,  Technical  Library)  (2  copies) 

Air  Force  Logistics  Center 
Robins  AFB,  GA  31098 

Attn:  WRALC/MMETE  (LT  W.  Shelton) 

Aerojet  ElectroSystems  Company 
A Div.  of  Aerojet-General  Corp. 

1100  W.  Hollyvale  Street 
Azusa,  CA  91702 

Attn:  A.  R.  Moorman 

AiResearch  Manufacturing  Co.  of  Arizona 
A Division  of  the  Garrett  Corp. 

P.0.  Box  5217 
Phoenix,  AZ  85010 

Attn:  G.  L.  Merrill 
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Armament  Systems,  Inc. 

712-F  North  Valley  Street 
Anaheim,  CA  92801 
Attn:  J.  Musch 

A.  T.  Kearney  and  Company,  Inc. 

100  South  Wacker  Drive 
Chicago,  IL  60606 
Attn:  R.  H.  Rose 

AVCO 

Lycoming  Division 
550  So.  Main  St. 

Stratford,  CT  06497 
Attn:  R.  Cuny 

Battelle  Memorial  Institute 
505  King  Ave. 

Columbus,  OH  43201 

Attn:  J.  H.  Brown,  Jr. 

The  BDM  Corp. 

2600  Yale  Blvd  SE. 

Albuquerque,  NM  87106 
Attn:  A.  J.  Holten 

The  BDM  Corp. 

1920  Aline  Ave. 

Vienna,  VA  22180 

Attn:  J.  W.  Milenski 

Beech  Aircraft  Corp. 

9709  E.  Central  Ave. 

Wichita,  KS  67201 

Attn:  Engineering  Library 

Bell  Helicopter  Textron 
Division  of  Textron  Inc. 

P.0.  Box  482 

Fort  Worth,  TX  76101 

Attn:  Security /Dept.  12,  J.  R.  Johnson 

The  Boeing  Aerospace  Company 
P.0.  Box  3999 
Seattle,  WA  98124 

Attn:  J.  G.  Avery,  M/S  41-37 
Attn:  R.  J.  Helzer,  M/S  47-28 
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The  Boeing  Company 
Vertol  Division 
Boeing  Center 
P.0.  Box  16848 
Philadelphia,  PA  19142 

Attn:  J.  E.  Gonsalves,  M/S  P32-19  (2  copies) 

Boeing  Wichita  Company 
3801  S.  Oliver  St. 

Wichita,  KS  67210 

Attn:  H.  E.  Corner,  M/S  K16-67 
Attn:  D.  Y.  Sink,  M/S  K16-14 

Booz. Allen  Applied  Research 
4330  East  West  Highway 
Bethesda,  MD  20014 
Attn:  W.  Djinis 

Calspan  Corp. 

P.0.  Box  235 
Buffalo,  NY  14221 

Attn:  Library  (V.  M.  Young) 

COMARCO  inc 
1417  N.  Norma 
Ridgecrest,  CA  93555 

Attn:  G.  Russell  (2  copies) 

Denver  Research  Institute 
University  of  Denver 
University  Park  Station 
2360  S.  Gaylord  Street 
Denver,  CO  80210 

Attn:  C.  R.  Hoggatt 
Attn:  R.  F.  Recht 

ETI  Effects  Technology,  Inc. 

A Subsidary  of  Flow  General,  Inc. 

5383  Hollister  Avenue 
Santa  Barbara,  CA  93111 

Attn:  Library  Acquisitions  (S.  Clow) 

Fairchild  Industries,  Inc. 

Fairchild  Republic  Division 
Conklin  Ave. 

Farmingdale,  L.I.,  NY  11735 
Attn:  J.  A.  Arrighi 
Attn:  G.  Mott 
Attn:  D.  C.  Watson 

Attn:  Engineering  Library  (G.  A.  Mauter) 
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Falcon  Research  and  Development  Co. 

2350  Alamo  Ave. , SE 
Albuquerque,  NM  87106 
Attn:  W.  L.  Baker 

i 

Falcon  Research  and  Development  Co. 

696  Fairmount  Ave. 

Baltimore,  MD  21204 
Attn:  J.  A.  Silva 

Firestone  Tire  & Rubber  Co. 

Firestone  Coated  Fabric  Co.  Division 
P.0.  Box  869 
Magnolia,  AR  71/53 

Attn:  D.  L.  Byerley 

Ford  Aerospace  and  Communications  Corp. 

Ford  Road,  P.0.  Box  A 
Newport  Beach,  CA  92663 
Attn:  Library 

General  Dynamics  Corp. 

Convair  Division 
P.0.  Box  80986 
San  Diego,  CA  92138 

Attn:  Research  Library,  MZ  40-6540  (U.  J.  Sweeney) 

General  Dynamics  Corp. 

Fort  Worth  Division 
Grants  Lane,  P.0.  Box  748 
Fort  Worth,  TX  76101 

Attn:  P.  R.  deTonnancour/G.  W.  Bowen 

General  Electric  Co. 

Aircraft  Engine  Group 
1000  Western  Ave. 

West  Lynn,  MA  01910 

Attn:  J.  M.  Wannemacher 

General  Electric  Co. 

Aircraft  Engine  Business  Group 
Evendale  Plant 
Mail  Drop  H-9 
Cincinnati,  OH  45215 

Attn:  AEG  Technical  Information  Center  (J.  J.  Brady) 
Attn:  B.  Alexander  (2  copies) 
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General  Research  Corporation 

P.0.  Box  6770 

5383  Hollister  Ave. 

Santa  Barbara,  CA  93111 
Attn:  J.  H.  Cunningham 
Attn:  R.  Rodman 

General  Research  Corporation 

SWL  Division,  Suite  700,  Park  Place 

7926  Jones  Branch  Dr. 

McLean,  VA  22101 
Attn:  T.  King 

Goodyear  Aerospace  Corp. 

Engineered  Fabrics  Division 
1210  Massillon  Rd. 

Akron,  OH  44315 

Attn:  Library,  D/152G  (R.  L.  Vittitoe/J.  R.  Wolfersberger)  (3  copies) 

Grumman  Aerospace  Corp. 

South  Oyster  Bay  Rd. 

Bethpage,  NY  11714 

Attn:  J.  P.  Archey  Jr.,  Dept.  662,  Mail  C42-05 
Attn:  R.  W.  Harvey,  Mail  C27-05 
Attn:  H.  L.  Henze,  B27-07 
Attn:  J.  Noack,  M/S  BIO-25 

Attn:  Technical  Information  Center,  Plant  35  LOl-35  (H.  B.  Smith) 

Hughes  Helicopters 
A Division  of  Summa  Corp. 

Centinela  Ave.  & Teal  St. 

Culver  City,  CA  90230 

Attn:  Library,  2/T2124  (D.  K.  Goss) 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Attn:  Technical  Information  Center,  DIMO  (P.  Okamoto) 

IIT  Research  Institute 
10  West  35  Street 
Chicago,  IL  60616 
Attn:  I.  Pincus 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  MD  20810 
Attn:  C.  F.  Meyer 
Attn:  V.  W.  Woodford 
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Kaman  Aerospace  Corporation 
Old  Windsor  Rd. 

Bloomfield,  CT  06002 
Attn:  H.E.  Showalter 

Lockheed-California  Co. 

A Division  of  Lockheed  Aircraft  Corp. 

2555  Hollywood  Way 
P.0.  Box  551 
Burbank,  CA  91520 

Attn:  Technological  Information  Center,  84-40  Unit  35,  Plant  A-l 
Attn:  G.  E.  Raymer,  D/75-84  Bldg.  63  A-l  (2  copies) 

Attn:  A.  D.  Jackmond,  Dept.  75-60,  Bldg.  170  B-l 

Lockheed-Georgia  Co. 

A Division  of  Lockheed  Aircraft  Corp. 

86  S.  Cobb  Drive 
Marietta,  GA  30063 

Attn:  W.  T.  Mikolowsky,  72-08  Zone  415 

Attn:  Sci-Tech  Info  Center,  72-34  Zone  26  (T.  J.  Kopkin) 

Lockheed  Missiles  & Space  Company,  Inc. 

A Subsidiary  of  Lockheed  Aircraft  Corp. 

P.0.  Box  504 
Sunnyvale,  CA  94088 

Attn:  G.  R.  Evans  (5501-572-5) 

Martin  Marietta  Corp. 

Orlando  Division 
P.0.  Box  5837 
Orlando,  FL  32855 

Attn:  Library  (M.  C-  Griffith,  MP-30) 

McDonnell  Douglas  Corp. 

Douglas  Aircraft  Company 
3855  Lakewood  Blvd. 

Long  Beach,  CA  90846 

Attn:  Technical  Library,  Cl-250/36-84  AUTO  14-78  (3  copies) 

McDonnell  Douglas  Corp. 

P.0.  Box  516 

St.  Louis,  MO  63166 

Attn:  R.  D.  Detrich,  Dept.  022 

New  Mexico  Institute  of  Mining  and  Technology 
Campus  Station 
Socorro,  NM  87801 
Attn:  Tera 


51 


JTCG/AS-78-V-004 


Northrop  Corp. 

Aircraft  Division 
3901  W.  Broadway 
Hawthorne,  CA  90250 

Attn:  J.  H.  Bach,  2130/83 
Attn:  V.  B.  Bertagna,  3451/81 
Attn:  H.  W.  Jones,  3360/82 
Attn:  J.  F.  Paris,  3628/83 

Northrop  Corp. 

Ventura  Division 

1515  Rancho  Conejo  Blvd. 

P.0.  Box  2500 
Newbury  Park,  CA  91320 
Attn:  M.  Raine 

The  Rand  Corp. 

1700  Main  St. 

Santa  Monica,  CA  90406 
Attn:  N.  U.  Crawford 

R&D  Associates 
P.O.  Box  9695 
Marina  Del  Rey,  CA  90291 
Attn:  B.  Jaeger 

Rockwell  International  Corp. 

Los  Angeles  Division 
International  Airport 
5701  W.  Imperial  Hwy 
Los  Angeles,  CA  90009 
Attn:  R.  Hurst,  MB56 
Attn:  W.  L.  Jackson 
Attn:  S.  C.  Mellin,  AD-25 
Attn:  R.  Moonan,  AB78  (2  copies) 

Attn:  T.  C.  Getten,  MB56 

Rockwell  International  Corp. 

4300  E.  Fifth  Ave. 

P.O.  Box  1259 
Columbus,  OH  43216 

Attn:  Technical  Information  Center  (D.  Z.  Cox)  (2  copies) 

Science  Applications,  Inc. 

200  Lomas  Blvd. , N.  W. , Suite  1020 
Albuquerque,  NM  87102 
Attn:  Library 


52 


JTCG/AS-78-V-004 


Science  Applications,  Inc. 

2361  Jefferson  Davis  Hwy. 

Arlington,  VA  22202 
Attn:  D.  A.  Venor 

Sikorsky  Aircraft  Division 
United  Technologies  Corp. 

North  Main  Street 
Stratford,  CT  06002 

Attn:  D.  P.  Bartz,  Chief  Survivability 

Southwest  Research  Institute 
P.0.  Drawer  28510 
San  Antonio,  TX  78284 
Attn:  Bessey,  Div.  02 
Attn:  P.  H.  Zabel,  Div.  02 

Stanford  Research  Institute 
333  Ravenswood  Ave. 

Menlo  Park,  CA  94025 
Attn:  J.  Golins 

System  Planning  Corporation 
1500  Wilson  Blvd. , Suite  1300 
Arlington,  VA  22209 
Attn:  J.  A.  Navarro 

Telcom  Systems,  Inc. 

A subsidiary  of  Telcom  Incorporated 
2300  South  9th  St. 

Suite  400 

Arlington,  VA  22204 
Attn:  C.  Gentzel 

Teledyne  Ryan  Aeronautical 
2701  Harbor  Dr. 

San  Diego,  CA  92112 

Attn:  Technical  Information  Services  (W.  E.  Ebner) 

Uniroyal  Inc. 

312  N.  Hill  Street 
Mishawaka,  IN  46544 
Attn:  J.  D.  Galloway 

United  Technologies  Corporation 
Government  Products  Division 
P.0.  Box  2691 

West  Palm  Beach,  FL  33402 

Attn:  P.  E.  Desrosiers,  Security  Officer 
Attn;  J.  Fyfe,  Mail  E-39 
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University  of  Dayton 
Research  Institute 
300  College  Ave. 

Dayton,  OH  45409 

Attn:  M.  Goldschmidt 

Vought  Corporation 
P.0.  Box  225907 
Dallas,  TX  75265 

Attn:  D.  M.  Reedy,  2-30100 

Williams  Research  Corp. 

2280  W.  Maple  Rd. 

Walled  Lake,  MI  48088 
Attn:  Library 
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